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Abstract 



£^ \ We discuss the possibility, that multiply charmed exotics are produced in 

t^j- ' ultrarelativistic heavy ion collisions. Production probabilities at RHIC are 

estimated and found to be large enough to allow for the possible detection 
of single and double charmed exotics. We demonstrate that charm matter is 
bound due to an attractive one-gluon exchange potential. Furthermore, we 



study the color magnetic and color electric potential separately, as well as the 
group structure of such exotics and estimate their masses. 



g \ I. INTRODUCTION 

S3 

Relativistic heavy ion collisions at very high energy will produce several charm quarks 
in one single event. Estimates using perturbative QCD predict about 10 produced charm- 
anticharm quark pairs at the Relativistic Ion Collider (RHIC) at Brookhaven in one single 
central collision of two led nuclei at an bombarding energy of \fs = 200 AGeV. This offers 
the unique opportunity to study the possible existence of multiple charmed objects. 

It is remarkable, that hadrons with two or more charm quarks, like the Q cc and the 
Q ccc have not be seen so far (for a summary see e.g. jl]). Certain exotic states with charm 
quarks have been discussed intensively in the literature. Multiquark states have been stud- 
ied in potential models using meson exchange • The tetraquark with two heavy quarks 
of mass rriQ and two light (u,d,s) quarks is bound in the limit of rriQ — > oo, but the sit- 
uation for the charm quark with m c « 1.5 GeV is less clear (see e.g. || and references 
therein). A recent investigation finds that the heavy Tetraquark system is bound when 
using pseudoscalar-meson exchange contributions ||. The pentaquark with one anticharm 
quark and four light quarks was proposed independently by Lipkin |7|] and by Gignoux and 
coworkers ||. This state is bound by an enhanced color magnetic binding energy. Hex- 
aquark and Heptaquark states with charm have been discussed also |J. There are proposals 
to search for the tetraquark and pentaquark in pp collisions (for an overview see e.g. []10| , pT| ) 
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at CHARM2000 and CERN COMPASS. Other exotics, especially for more than two charm 
quarks, have not been considered so far. 

Here, we investigate the structure of multiply charm objects (charmlets), as well as their 
masses and stability. First, estimates are given for the production rates of these objects at 
RHIC using a simple coalescence approach. Then, the properties of charm bulk matter are 
discussed. A modified bag model is introduced in section three which is able to describe 
hadrons with and with out charm in the same approach The color electric and magnetic 
interactions for charmlets are investigated in section four. Finally, we give a summary and 
an outlook. 



II. PRODUCTION OF CHARMED EXOTICS AT RHIC 

The first and most basic question we have to address is, if there is any chance to de- 
tect multi-quark exotics in heavy ion collisions. We will try to answer this question in this 
section in the case of Au + Au collisions at RHIC. First we estimate very crudely the pro- 
duction probability of multi-charmed hadrons and then compare this production probability 
to detailed calculations within the cascade model HIJING Jl2| using a coalescence code as 



afterburner. An early estimate of heavy baryon production at RHIC and LHC within a 
simple statistical model can be found in ||13|| . 

The first scale we have to determine is the total number of events N events at the RHIC 
collider in an average year. Every production rate or probability P has to be compared to 
this number. Only in cases for which N events P > 1 we have a fair chance to measure the 
exotic. Given the luminosity at RHIC ~ 2 10 26 cm -2 s -1 and the geometric cross section of 
gold nuclei ~ 3 barn, we find that there can be up to 2 10 10 heavy ion collision events per 
year. This number has to be reduced by several factors. First, the accelerator will run only 
part of the year. In addition, the detectors will only cover a small fraction of the full solid 
angle and one has to face all kind of inefficiencies. Finally, there is the branching ratio of 
the exotic state to measurable (charged) final states. 

Guided by the proposal for the experiments CHARM2000 and CERN COMPASS de- 



signed to measure charmed exotics in pp collisions |T0| , pT| , one might be able to recover as 
little as 0.1% of the actual data. Combining these reduction factor with the integrated lu- 
minosity for a month or so, we obtain a very conservative estimate for the expected number 
of events at RHIC per year N events « 10 6 . We can expect, that an experimental observable 
occuring with a likelihood larger than 10~ 6 can be detected. 

A simple method to estimate the production probability of charmed hadrons is the coa- 
lescence model [|14|| . In its basic version this model implies, that the production probability 
of a multi-quark state should be proportional to the product of the ratios of the number 
densities of its constituents normalized to the total density of the surrounding matter. In 
hot and dense hadronic matter charmed quarks are mainly formed via pair creation. The 
cross section for this pair creation process for proton nucleus collisions has been measured 



T5| and is in the order a C c ~ 0.04 mb for a lab energy of Ei a b = 800 GeV. Extrapolated 



to the bombarding energy available at RHIC of ^/s = 200 GeV one gets a c5 ~ 300 fib |T6| 



This amounts to about 10 cc pairs in one single central collision of gold nuclei at RHIC! 
The charm pair production cross section has to be compared to the total nucleon nucleon 
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cross section of <jnn ~ 40 mb. As a first guess, we find that an exotic containing C charmed 
quarks is produced with a probability proportional to (ucc/^nn) ~ 1CT 2C . 

This guess neglects the light quark content of the exotic. Actually, in equilibrated 
hadronic matter the charmed quark-anti-quark pairs would annihilate again and we would 
end up without any charmed hadrons. Only the surrounding light quarks can combine 
with the charmed quarks to form D-mesons, £ c -baryons etc. In this way we maintain 
the charmed quarks, prevent their annihilation and enhance the probability that these type 
of singly charmed objects combine with itself or other hadrons to multi-charm and multi- 
baryon number carrying exotics. High light quark density is therefore essential to assure the 
production of these kind of exotics. The central rapidity region should be the best area to 
find them experimentally. 

An estimate of both, the light quark and baryon number density, is the number of 
nucleons iVjv of the collision normalized by the number of pions N n expected for collisions at 
RHIC, i.e. several thousand. The probability to find an charmed exotic with baryon number 
A is therefore proportional to the ratio (Njy/N n ) A . At central rapidity this might be around 
10~ 2A . The suppression factor for forming a deuteron is then about 10~ 4 . 

Putting these results together, we find that the production probability of a charmed 
exotic with charm number C and baryon number A is 

P(C,A) = K ( (^) A « io^- 2 A (1) 

From this estimate we gain a first feeling what possibly could be detected at RHIC. Eq. (]!]) 
gives one deuteron in 10 events and 10 D-mesons per event. Combined with the supposed 
sensitivity of 10~ 6 the inequality 

2A + 2C<9 (2) 

for the detectability of multiply charmed exotics emerges. A single charmed exotic carrying 
baryon number three, a double charmed exotic with baryon number two, or a triple charmed 
exotic with baryon number one (the fi ccc , the charmed analog to seems to appear in 
enough events to be visible in the detectors. 

Let's solidify this result with a simulation. For a first estimate we used the momentum 
distribution of charm quarks from the event generator HIJING |TJj as input for a simple 
coalescence model. Charm quarks sitting in within a certain radius of momentum p are 
supposed to form a multiply charmed exotic. We assume, that the main penalty comes from 
the coalescence of charm quarks and neglect the penalty of finding light quarks within the 
same area as they are abundantly produced at such high energy. Recent investigations of 
charm quark production at RHIC indicate a sizable energy loss of charm quarks during the 
expansion stage especially in transverse momentum. A so called quench incorporates the 
energy loss of the charm quarks by rescattering with the surrounding particles after being 
produced in a heavy-ion collision. This makes the phase space distribution of the charm 
quarks narrower than in the unquenched case (for a more detailed discussion about this point 
see [12]). Fig. |3] shows the number of produced charmed exotics for various charm numbers 
as a function of the coalescence parameter p . Here we choose the case with shadowing and 
a quench, the case for no shadowing without a quench is shown in Fig. [|. In both cases, the 
production of double and triple charmed objects seems to be detectable at RHIC, i.e. the 
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production rate is above 10~ 6 per event. A quench enhances the formation of charmed matter 
as the heavy quarks are closer together in phase space. We find indeed, that for the quenched 
case, the production rate is higher and less sensitive to the coalescence parameter p as the 
charm quarks. Another scenario was recently discussed in [|17| which will also enhance the 
production rates of charmed matter droplets. Within a cascade hydrodynamical model, they 
find that charmed particles feel drag and diffusion forces within quark-gluon plasma droplets 
so that the charm quarks resides inside the droplet. This mechanism resembles very much 
the separation mechanism for strange quarks: the quark-gluon plasma droplet gets enriched 
with strangeness as it is energetically more favorable to keep the strange quark inside the 
droplet . Here we point out that the distillation mechanism proposed for the formation of 
strangelets Jl9[ will also work for charmed matter, as the D-meson is much lighter than the 
lightest baryon with charm A c . This effect will furthermore enhance the charm fraction of 
the droplet during cooling down by emitting D-mesons and carrying away anticharm quarks, 
therefore enriching the plasma droplet with charm. Nevertheless, this implies a baryon-rich 
regime or strong fluctuations of baryon number which one might expect at RHIC energy 
PDfl . Therefore, charmed exotics of \C\ > 3 might well be produced in measurable amounts 
if a quark-gluon plasma forms. Hence, we will investigate in the following the properties of 
charmed objects for \C\ > 3, too. 



III. CHARM MATTER 

First, let us discuss the properties of charm matter in bulk. Strange matter gets enhanced 
stability compared to two flavor matter by introducing a new degree of freedom and a 
lowering of the Fermi energy pl"| . This is possible, as the strange quark mass is still light 
enough. For the charm quark with a mass of the order of 1.5 GeV, the Fermi energy can not 
be lowered at reasonable density. Calculations for hybrid stars with charm quarks conclude 
that this quark is probably too heavy to appear in the core of a neutron star [p2| |. On the 



other hand, the potential generated by an effective one-gluon exchange gets attractive and 
gives an enhanced stability for heavy quarks. This potential is repulsive for massless quarks 
and is evident from the energy density for massless quarks which reads 



7T 



Here a c is the strong coupling constant and Hq is the chemical potential of the up and down 
quarks. Note that for values of a c > tt/8 ~ 0.39 the correction from one-gluon exchange 
is bigger than the zeroth order. Usually, bag model fits to hadron spectra yield a c > 7r/8 
23| , |24|1 indicating a highly non-perturbative potential. This is a problem for bulk matter, as 



the pressure gets negative because 

Therefore, the parameters extracted from a fit to hadron spectra can not be used for bulk 
matter calculation. This indicates that the bag model is an effective model and that the 
parameters change when going to higher density and/or strangeness numbers. 
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The situation changes for a finite quark mass. The energy density for quarks with a 
finite mass is given by [25 



8tt 2 



fii + X 16a c J / 2 fii+X\ 2 / 2 2 \ 2 



0» ( 2/Xi -m i )-m i log + — — < 3 [mx - m x log = j -2[fjq-mi 

rrii Sir V rrij 



(5) 



with x = \J [if — mf. Here the interaction gets attractive for large enough quark mass m*. 
This can be seen in Fig. 0, where the total energy per baryon number (minus the mass) is 
plotted versus the charm fraction f c = \C\/A for increasing values of the strong coupling 
constant a c . For a c = 0.3, charm matter is bound by 55 MeV at f c = 1.7. Hence, the 
charm quark is heavy enough to feel an attractive interaction from the effective one-gluon 
exchange. For comparison, the case for strange matter indicated by m s is also shown, then 
as a function of the strangeness fraction f s = \S\/A. Only the case a c = can be seen on 
the plot. Adding the one-gluon exchange contribution results in a repulsive term, so that 
the curves for the strange quarks are shifted up and are lying outside the range displayed in 
the figure. 

In general, charm matter consists of strange quarks also. We plot the areas, where charm 
matter is bound in a contour plot as a function of the charm and strangeness fraction in 
Fig. |2| for different values of the strong coupling constant a c . Note that charm matter can 
not exist for f s + f c > 3 by definition (one can not have more than three quarks per baryon 
number). The calculation was done for a bag parameter of B l / A = 235 MeV taken from a fit 



to charmonium levels jZfJ which is in accordance with estimates from QCD sum rules 
Also shown is the case for the parameters taken from a fit to the hadron spectra (denoted by 
fit) which we will discuss in the following section. This is the area in the lower right corner 
of the figure with a much lower bag parameter of B 1 / 4 = 168 MeV. Note that the pressure 
is negative for the massless quarks in this case as a c = 0.476 > 7r/8! 

Charm matter is bound when the total energy per baryon is lower than the sum of the 
corresponding hadron masses. The masses of the hadrons with C > 2 are not known and 
taken from the bag model calculation (see Table |J). Strange matter along the line for 
f c = is unstable for this bag parameter [I9| except for the parameters taken from the 



fit. Nevertheless, charm matter demonstrates to be stable for a wide range of charm and 
strangeness fraction of f c = 1 — 3 and f s = — 1.5 even for this high bag parameter. Hence, 
charm matter might be stable while pure strange matter is not. The next step is to calculate 
the masses of finite charm matter for low baryon numbers. This will be done in the following 
sections after introducing the modified bag model for multiply charm quark bags. 



IV. MODIFIED BAG MODEL 

At present our technical skills in solving QCD do not suffice to accurately predict bound- 
state energies of multi-quark states. Instead we have to rely on simple model calculations 
and estimates. A very simple, controllable and relatively reliable model for our purposes is 
the bag model p3 |. We think the results of this model give us a qualitative first impression 



of the stability of multi-quark states. More detailed calculations can be performed both 
through potential interactions and on the lattice. 
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We modify the commonly used version of the bag model in two essential ways. First, 
we evaluate the bag energy by setting the strong coupling a c to zero effectively 
neglecting all color electric and magnetic contributions. Both contributions are evaluated 
after minimization of the bag energy, then added to the bag mass, and therefore treated 
as true lowest order corrections. In this way we take only the gluon exchange contribution 
into account and our bag is locally color charged, while globally color neutral. We postulate 
that all other color electric and color magnetic corrections can be accounted for by adjusting 
the parameters of the model. And secondly, we treat all quarks on the same footing, i.e. 
we do not distinguish between heavy and light quarks, rather demand flavor symmetry 
which is broken by the explicit masses of the quark species. This simplifies our approach. 
A very heavy quark, like the bottom is thus not treated non-relativistically and is not 
sharply localized at the center of the bag. Instead it is distributed, with the maximum 
of the distribution being closer to the center of the bag than the one of the light quarks. 
Both of these modifications are ad hoc and can only be justified by the success of the 
phenomenological model in describing the experimentally measured hadronic masses. 

We are interested in a first crude estimate of the binding energy of a multi-quark state 
with n quarks. All quarks are confined inside a spherical bag of radius R. The energy or 
mass of such a state is to zeroth order in the strong coupling a c given by 

n Z 4 

E n = J2^-^ + -nBR 3 , (6) 
i=x K 6 

where u>i is the eigenenergy of the z-th quark confined within the bag, Z is the Casimir 



energy of the bag [^5 absorbing all effects ~ l/R and B is the bag constant. Minimization 
of the bag energy determines the bag radius R. 

The confined i-th quark has to satisfy both the free Dirac equation and the linear bound- 
ary condition at the bag radius R which leads to the constraint for the single particle energy 
Xi of a quark 

tanxj = — -, (7) 

1 - R (uji + mi) K ' 

where rrii is the mass and u>i the total energy. In Fig. |5] we plot the probability density p of 
a quark of flavor q(u, d), s, c, b for a bag of radius R = 1 fm as a function of the scaled radius 
r/R. The probability density is proportional to the sum of the square of the spherical Bessel 
functions 

p(r) = ^(r,t) t 7o ^(r,t) 
N 2 

= — ({uo + m)] Q (x i r/R) 2 + (uj-m)] l {x i r/RY) . (8) 

The heavier the quark, the more localized is the density at r = and the smaller are the 
values for the spherical Bessel Functions. 

The minimization of the bag energy in eq. (|6|) fixes the bag radius R. We use this radius 
to determine the color electric and color magnetic corrections to order a c . Both corrections 



were discussed by many authors, e.g. p3| , p4] , |28| and we quote only the results. The color 
electric correction Ve is 
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i>j a JO u 

where the A" are the color matrices of particle i and the density pi is given by 

, cjj [scu — sin 2 (a;M)/a;M] — m* [sin(xu) cos(xu) — sm 2 (xu)/xu] 
u)i [x — sin 2 (a;)/x] — mi [sin(x) cos(x) — sin 2 (x)/x] 

Note, that the summation is performed over all quark pairs in the bag (i > j). In the 
original MIT bag model, the sum runs over all quark combinations |[23|| . Then, the 

color electric potential for the charmonium state will vanish 

= £ £(A? • >< E cc = C 3 [l] x E cc = (11) 

i,3 a 

as it is a color singlet state for which the Casimir operator of color SU(3) Cs[l] = 0. 
Therefore, the original MIT bag model can not describe the mass of J/^/, as it is mainly 
bound by the Coulomb-like color electric potential. In our case, the sum runs only over 
unequal quark pairs (i > j), therefore, the color electric potential turns out to be 

Vco = ££(A? • A") x E cc = f~C 3 [l] - ^) x E cc = - 1 4e cc (12) 



i>j 



a 



i.e. nonvanishing and strongly attractive, as needed to describe the masses of J/^f and r) c 
correctly. This procedure allows now for a local color charge (there is an intrinsic color 
charge inside the bag), but the overall color charge still vanishes as the sum over all quarks 
in the bag gives 

£££(^-\ a ) = £^3=cr a ' = o . (is) 

k i,j & k 

Here Cf is the Casimir operator of the quark flavor k = q(u, d), s, c. Hence, light quarks and 
heavy quarks can have a color charge and do not need to be color neutral separately, but 
of course the sum of all color charges gives an overall color neutral object. This new color 
degree of freedom will be important for the enhanced stability of multiquark states which 
we will discuss later in more detail. 
The color magnetic correction Vm is 



Vw = -££fa\ )-(°i*i)Mfc 

i>j a 

3a c fj,(rrii,R) n(rrij,R) ( f R dr 

i>j 



£ £(^A°) ■ (o- 3 - A?) R g (1 + 2 I ^M™^)/iK,r)j , (14) 



where the <7j are the spin matrices and the magnetic moment of the quark eigenmode 
given by 
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1 AruJi + Irrrii - 3 

uAmi.r) = ; r . (15) 

p v ' ' 6 2u; i (ra; i -l)+mi V 7 

The summation runs again over all quark pairs in the bag. 

In Fig. [3] we plot the electric interaction and the magnetic interaction between 
two quarks with flavor i,j = q(u,d),s,c as a function of the bag radius R. We see that 
the largest contribution to the total electric interaction energy within the bag is given by 
the interaction between two heavy quarks. In contrast, the largest contribution to the total 
magnetic interaction energy within the bag is given by the interaction between two light 
quarks. A multi-quark bag becomes therefore more stable the more heavy charmed quarks 
we add, due to the increase in the electric interaction energy, and it becomes more stable 
the more flavor antisymmetric we arrange the light quarks (u,d,s). This property will allow 
us to simplify the study of the group structure of multi-quark bags in the next section. 

As a test of our modified bag model we calculate first the measured masses of the light 
hadrons built up from u, d and s quarks. First, we use some of the known masses to fix the 
parameters of the model, 5 in the case of the light hadrons. The mass of the lightest quarks 
is set to zero m q = 0, where q = u,d. The strong coupling constant a c = 0.4764 is fixed 
from the color magnetic mass splitting of the proton p and the delta A. These two particles 
have identical color electric interactions . The masses of the proton and the mass of the 
degenerate p, cj-system are used to fix the bag constant B 1 ^ = 0.1684 GeV and Zq = 1.585. 
Finally, we use the Q~ hyperon to fit the strange quark mass m s = 0.342 GeV. The various 
color and spin assignments of the color electric and magnetic energies in @ and ( |i~4"D are 



given in [p3" |. 

In Table | we list our results for the light hadrons. We reach satisfactory agreement for 
the heavier hadrons, but have problems with the lightest mesons like the ir, rj system and light 
strange mesons like the K and the K*. We attribute this to their role as Goldstone bosons 
of chiral symmetry. The slightly larger error of the <fi might be due to its non-spherical 
shape. 

In the next step we include the two heavier quarks and check if our model continues to 
perform reasonably well. The charm mass m c = 1.788 GeV and the bottom mass mj, = 5.296 
GeV are fitted to the mass of the A c and the B s respectively. With these quark masses fixed 
we can predict all other hadrons containing u, d, s, c or b quarks. We find in Tables [Tl] and |IT1| 
good agreement between the experimental measured masses and our masses. In most cases 
we have a relative deviation from the measured masses of less than 2 — 3%. The predicted 
masses of Q cc and Q ccc are a little bit higher than estimates from potential models . With 
this in mind we can proceed in the next section to analyze the color and spin assignments 
for multi-charmed exotics and use then our modified bag model to predict their masses. 



V. POSSIBLE CANDIDATES AND THEIR MASSES 

A. Estimate of the color magnetic potential 

As seen in the previous section, quark bags can be bound by either the color magnetic 
force or the color electric force. The former one is dominant for light quarks, while the latter 
one is most pronounced for the heavy quarks c and b (see Figs. |] and 0). For charmlets, one 
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can show that both forces are attractive! Let us start with the color magnetic term. Here 



we follow the discussion outlined by Jaffe |31] for strangelets but with an unprecedented 



conclusion about charmlets. In this subsection, we consider up, down and strange quarks as 
light quarks and discuss effects due to the finite strange quark mass later. 

Three light quarks can be combined in SU(6) color-spin to a [56]-plet, two [70]-plets and 
a [20]-plet. A decomposition into SU (3) color x SU(2) spin shows that there is a color singlet 
for the [70]-plet (the baryon octet) and for the [20]-plet (the baryon decuplet), but not for 
the [56]-plet. The color magnetic term is proportional to 



V M (d 6 , d 3 , S, N q ) = £ ■ A") (aj ■ A?) Mij = (sN q - \c & [n] + ]-C 3 + ^S(S + 1)) M tJ 

i>j a \ I I 6 J 

(16) 



where N q is the number of light quarks, C% is the Casimir operator of SU(6) which depends 
on the representation [//], C 3 is the Casimir operator of SU (3) color, and S denotes the total 
spin. is the matrix element defining the strength of the interactions. The potential 
depends on the dimensions of color-spin SU(6) and color SU(3) following the notation of 
||. From eq. © one sees, that the representation with the highest Casimir operator Cq 



will be mostly bound. Generally, the representation which is mostly symmetric in color-spin 
will have this feature [Q. It turns out that the [56]-plet is the most symmetric and will 
be more bound than baryons. But this one is color charged and does not exist in nature 
then. Now we include heavy quarks, like the charm quark, which will have a negligible 
contribution to the color magnetic term due to its heavy mass. Then one can make the 
[56]-plet of light quarks color neutral by adding e.g three charm quarks (or antiquarks) and 
gain this additional binding energy. This is of course in line with the discussion about 
the hyperfine interaction for pentaquark states introduced by Lipkin |7]] but from a slightly 
different point of view. One can estimate the binding energy by assuming that the matrix 
element Mjj is the same for all light quarks (i.e. setting m s = m U( i) and independent of the 
radius of the bag. Then one gets for the binding energy 

B(N q = 3) = V M (56, 8, 1/2, 3) - V M (70, 1, 1/2, 3) = -UM + 8M 

= -6M=-6x( M *- M ^-110MeV 

16 



The necessary values for the Casimir operators can be found in []32j| . We note that the [56]- 
plet is a singlet in flavor, as overall Fermi statistics require that the most symmetric one in 
SU(6) color-spin must be most antisymmetric in SU(3) flavor. Hence, the proposed exotic 
states would have a quark content of {udsccc} (the H ccc dibaryon) and {udsccc} (similar to 
a bound state of two D- mesons and one D s meson). 

The same line of argumentation can be repeated for four light quarks, which are in 
a [210]-plet and are color neutral by adding an anticharm quark (this is the well known 
pentaquark P cs 0) or two charm quarks (the H cc dibaryon with {uudscc}, {uddscc}) or 
={udsscc}). The binding energy is then 

B(N q = 4) = V M (2W, 3, 0, 4) - V M (70, 1, 1/2, 3) = -16M + 8M 
= -8M = — (M A - M N ) w -150 MeV 
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The H cc has been also discussed in || . For five light quarks, one gets 

B(N q = 5) = V M (420, 3, 1/2, 5) - V M (70, 1, 1/2, 3) - V M (21, 3, 0, 2) = -16M + 16M = 

using C*6[420] = 358/3 which we computed from the Tables in [33| and the formula given 
by Jaffe |35[| . This does not support the possibility for the bound exotics H c . For the 
heptaquarks {uuddscc}, {uudsscc}, {uddsscc} the situation is different as 

B\N q = 5) = 14/(420, 3, 1/2, 5) - V M (70, 1, 1/2, 3) = -16M + 8M = -8M « -150 MeV 

as the two charm antiquarks combine to two D-mesons without any gain in the color magnetic 
energy. The case for N q = 6 gives 

B(N q = 6) = Vm(490, 1, 0, 6) - 2 x V M (70, 1, 1/2, 3) = -24M + 16M = -8M « -150 MeV 

The line of argumentation stops here as it is just the famous case for the H-dibaryon |3B| 
which is color neutral by itself. Hence, the exotics {uuddssccc} and {uuddssccc} only exist, if 
the H-dibaryon is bound and/or the interaction between heavy and light quarks is attractive 
enough. 



B. Estimate of the color electric potential 

The discussion of the color electric term is simpler, as it involves only color SU(3). 
Between the charm quarks one finds the color electric potential 

V E (d 3 , N c ) = £ £(A? • A?) E l3 = (\ct - \n) E cc (17) 

i>j a \Z 6 / 

where E cc is the matrix element for the color electric interaction. This is the Coulomb-like 
potential used in e.g. charmonium potential models. This potential increases linearly with 
the number of heavy quarks hence giving an enhanced stability for multiply charmed objects! 
In addition, we calculate also the color electric potential between the charm quark and the 
light quarks: 

/N q +N c N q Nc \ 

V^(d 3 , N q , N c )=[ £ ■ \j) - £ £(A? ■ A") - £ £(A? ■ A?) E l3 = -C c 3 E qc (18) 

\ i>j a i>j a i>j a J 

as the contribution is quite sizable (see Fig. |^). Interestingly, this interaction does not 
depend on the number of light and heavy quarks. Its strength is determined by the internal 
color structure and it is always attractive. Note that the quark bag is still globally color 
neutral but locally color charged. 

The matrix element increases with increasing mass as evident in Fig. |5] so that the 
tetraquark with two heavy quarks Q and two light antiquarks is bound in the limit of 
m,Q — > oo. The case for the tetraquark with charm (ccqq) is still under debate but the 
bottom quark seems to be heavy enough to support a bound state 0. 
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C. Candidates for the flavor SU(3) symmetric case 



With the preceding discussions about the main color magnetic and electric potential 
terms, one can now look for possible candidates in the flavor SU(3) symmetric case, ignoring 
a finite mass of the strange quark mass. 

We calculate the coefficients in front of the matrix elements for the magnetic interaction 
between the light quarks, for the electric interaction between the heavy quarks and the 
one between heavy and light quarks. Bags up to 6 light quarks and 6 heavy quarks are 
considered. The result is shown in Table |V| for charm quarks and in Table |V| for the case 
of anticharm quarks. The matrix elements are compared to the possible decay products on 
the right side of the tables. We consider also the two possible color representations. The 
color-spin multiplets can be decomposed into SU(3) co i or x SU{2) spin as 



N q = 2 
N q = A 
N q = 5 
N q = 6 



[21] D (3,0) + (6,1) 
[210] D (3,0) + (6,1) + ... 
[420] D (3, 1/2) + (6, 1/2) + 
[490] D (1,0) + (8,1) + ... 



For N q = 3, we take the [56]-plet for the color octet and the [70]-plet for the color singlet 
representation, respectively. 

One notes, that for certain candidates, like {3q3c}, the coefficients for the potentials for 
the candidate and its decay product are just the same. Hence, these candidates are not likely 
to be bound. In four cases, only the color magnetic term is larger for the candidate than for 
the decay products ({cc4q},{6c6q}, the pentaquark {c4q}, and {6c6q}). Other candidates 
get additional binding energy from the magnetic term and one of the color electric terms 
depending on the color charge representation. For the cases of different color states for the 
charm quarks, one observes that the one with the higher dimension in color SU(3) gives less 
binding for the magnetic term and the electric term E cc but more binding for the electric 
term between the heavy and the light quarks E qc . As the latter term is the most dominant 
contribution to the overall binding energy (see Tables [TV] and [VJ) , it might well be that 
multiply charmed objects are bound in a higher representation of color SU(3)! 

To clarify this, we estimate the matrix element from Figs. |5] and ^ to M = 18 MeV, 
E qc = 50 MeV and E cc = 70 MeV and calculate the binding energy. Table [VI] shows 



the result. Indeed, for some cases the higher color charge representation turns out to be 
more stable: for {3c6q}, {cc2q}, {cc5q}, and {cc§q}. Accidentally, it comes out to be the 
same for two cases ({cc4q} and {4c4q}). Note, that we predict on this basis, that the 
tetraquark {cc2q} will be more bound, if the charm quark is sitting in the higher color 
SU(3) representation, i.e. in the color sexplet instead of the usually assumed color triplet. 
The tetraquark is then bound by the electric interaction between the heavy and light quarks 
and not by the Coulomb-like potential between the heavy quarks! The binding energy is 
quite sizable but will be diminished by flavor SU(3) breaking effects and finite size effects 
due to the different radii of the candidate and the decay products. 
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D. Flavor breaking and finite size effects 



For the flavor breaking case, one has to calculate the coefficient for every quark-quark 
combination. In the following, we take into account the finite strange quark mass and 
consider three different quarks species: q(u,d), s,c. 

The color electric contribution between the same quark species in the color state defined 
by C3 is given by 

N k I o 

e fcfc = EEW ■ A?) = -C7* - |iV fc (19) 

i>j a 

where the sum runs only over the quark species k. Hence, for the charmonium one gets 
e cc = —16/3, but for a baryon with two charm quarks it is e cc = C 3 [3]/2 — 16/3 = —8/3, 
a factor of two lower, often encountered in potential models. The case of three different 
quark species can be found by using the definition of the quadratic Casimir operator for 
color SU(3) 

N k +N t N k N t N k +N t 

E = + -a?) + 2 e E(\° ■ a-) = c*+* = c 3 m 

(20) 

where the Casimir operator of color SU(3) for the two quark species is just the one of the 
remaining quark species m due to color neutrality. This results in 



N k +Ni 

e k i = 

i>3 



k+Nl n 

E E^A°) = ^3 m -C 3 fc -Cl) (21) 



The same procedure can be done for the color magnetic interaction, here for SU(6). Defining 

««= E ( 22 ) 

i>j a 



gives 

a kk = 8N k - l -C k M + \ct + p(S + 1) (23) 
for equal quarks, otherwise 

oh = -\ (ct l - cl -c l 6 ) + \ (c 3 m - c$ - C l 3 ) + i {S{S + 1) - s k {s k + 1) - Sl ( Sl + 1)) 

(24) 

where S denotes the total spin of the subsystem of the quark species k and I. From Fig. [7] 
one sees, that one can neglect the color magnetic contribution for the charm quark. Then 
one can calculate the color magnetic interaction for the light quarks for the different possible 
cases of spin and color by taking care of the overall Fermi statistics. For example, for three 
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light quarks {uds} in a color octet state one gets Cf = CJ = 16/3, CJ = 12, a 



99 



a qs = — 6 or Cf = 40/3, CJ = 16/3, CJ = 12, a 99 = —4/3, a gs = —38/3, respectively. Here q 
stands for up- and down- quarks. For the color singlet case, only one combination is possible 
with Cf = CJ = 16/3, CJ = 0, a qq = —8 which is just the case for the A hyperon. With the 
coefficients determined above, we calculated the masses of all possible charmlet candidates 
within the modified bag model up to 6 light and 6 charm quarks. The result is shown 
in Table |V11| for charm quarks and in Table |V111| for anticharm quarks. Only in two case, 



indicated by a star, the threshold mass for the strong decay is known experimentally. All the 
other threshold values are estimated within the same bag calculation. As it is evident from 
the tables, none of the candidates has a mass lower than the (estimated) threshold mass, 
hence they are all unstable. The bag calculation might be accurate within 3% guided by our 
fit to the hadron spectra. There are four candidates which have a mass within this range of 
3% above the threshold mass, three with baryon number two: {cssudq} with charge +1 or 
0, {ccsudq} with charge +2 or +1, {ccssud} with charge +1, and one with baryon number 
three: {cccssuudd} with charge +2. None of the anticharm candidates comes so close to 
the threshold. The reason is that the decay products of the latter ones are D-mesons, which 
have a bigger binding energy than A c baryons due to their smaller radius. 

The dependence of the color magnetic and electric potential is quite sensitive to the radius 
of the bag. A larger bag feels therefore a smaller attractive potential than a smaller one. 
Secondly, the finite mass of the strange quark also decreases the color magnetic interaction. 
These effects can then result in an unbound candidate, which was bound in the flavor 
symmetric case. Nevertheless, there is another more important term, which determines 
whether or not a multi quark bag is bound. This is the so called Casimir energy which 



is proportional to —Z /R as also discussed in [p8| . As the decay products will involve at 
least two bags, this attractive term alone will raise the multi quark bag considerably. In the 
case of equal radii, say 1 fm, one gets an additional shift of Zq/R « 300 MeV! In reality, 
this shift will be even higher, as the radius increases proportional to the number of quarks 
sitting in the bag. Note, that this term is purely phenomenological and absorbs all effects 
proportional to 1/R for the hadron spectra. This might change considerably for multi quark 
states, but can not be determined within this simple bag model. All the charmed candidates 
below a mass of 8 GeV are lying within 300 MeV above the threshold. For the anticharm 
candidates, only the pentaquarks and the tetraquarks fulfill this constraint. 

Furthermore, the bag model assumes a spherical symmetric shape and an uniform radius 
for all quark species. This assumes, that the charm quarks are sitting in the core surrounded 
by the light quarks. As demonstrated in the previous sections, the color electric term between 
the light and heavy quarks might be the dominant one. This will cause the charmlet to form 
a highly nonspherical configuration, as the light quarks are located around each charm quark, 
while the charm quarks, feeling a repulsive potential between each other, are separated from 
each other. This structure has to be studied in a fully three-dimensional model which is 
beyond the scope of the present investigation. 

VI. SUMMARY 

We demonstrate that charm matter can be produced at the coming RHIC collider in 
reasonable amounts. We estimate the production rates of these exotic new form of matter 
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at RHIC in a simple coalescence approach and find, that about 10-100 multiple charmed 
exotics can be measured per year. Energy loss of the charm quarks and hydrodynamical 
effects can enhance these production rates considerably. 

We investigate the properties of charmed matter and finite multiple charmed objects. 
The one-gluon exchange in bulk matter gets attractive for a sufficiently high mass of the 
charm quark, much more than for the lighter strange quark. Charm matter is bound with 
respect to strong decay to hadrons within an area of charm fraction of f c = 1 — 3 and a 
strangeness fraction / s = — 1.5 even if strange matter is unstable. 

We introduce a modified bag model where the parameters are fitted to hadron masses. 
The model can describe the masses of hadrons with and without charm on the level of a few 
percent. The predictions of the masses for multiply charm hadrons are a little bit higher 
than the ones from potential models. 

We discuss then the color magnetic and color electric potentials for multi quark states 
with charm. For finite systems, we show that in the limit of SU(3) flavor symmetry the 
color magnetic interaction can be more attractive than for the corresponding hadron with 
a similar quark content. Systems with 3 and 4 light (u,d,s) quarks have an enhanced color 
magnetic term compared to ordinary hadrons when adding heavy quarks. The constraint of 
color neutrality can be released for the light quarks as the heavy charm quarks neutralizes 
the overall quark bag. This allows for an additional freedom for the spin and color structure 
of the light quarks, therefore enhancing the binding energy. In addition, the color electric 
interaction between the charm quarks is also attractive and increases with the number of 
charm quarks. It turns out that the dominant contribution to the binding of multi charm 
objects comes from the color electric term between heavy and light quarks which depends 
solely on the internal color structure. These effects result in an enhanced stability of multiply 
charmed exotics. 

For flavor symmetry and ignoring finite size effects, several charmlet candidates have 
binding energy of more than 100 MeV. A full calculation within the modified bag model 
does not give a bound candidate. Also the pentaquark and the tetraquark are not bound 
in our approach. For some other cases the masses are lying just a few percent (within 
our range of uncertainty, i.e. less than 3%) above threshold. These candidates are: the 
hexaquarks {cssuud}, {cssudd} with charge +1 and 0, {ccsuud}, {ccsuud} with charge +2 
or +1, {ccssud} with charge +1, and one candidate for C = 3 {cccssuudd} with charge 
+2. These candidates might be bound, if one replaces the charm quarks with the heavier 
bottom quarks, as the attractive color electric forces will increase for heavier quarks. Note, 
that the Tetraquark is bound in the limit of mq — > oo. We checked this for the above 
candidates but found that the candidates with bottom quarks are still lying above the 
estimated threshold. The main reason for the sudden instability of the candidates in the 
full bag model is, that the effective energy correction —Zq/R shifts multi quark states up 
by at least +300 MeV. Another reason is that non-spherical solutions might exist which are 
lying lower in energy. This picture is supported by the strongly attractive potential between 
light and heavy quarks which will cause the charm quarks to be separated and not sitting 
in the center of the bag as assumed in our calculation. This can only be addressed in full 
three-dimensional calculation. Further investigations along this line using other models, like 
potential models or lattice calculations, are certainly needed to clarify this question. In 
this context, it is interesting to note that meson exchange contributions are important for 
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binding multiquark systems as pointed out by several authors 

As the charmed hadrons live on a timescale around r »s lOOyU a vertex detector is needed 
for their detection which must be sitting extremely close to the target. Micro-strip detectors 
are used for the CHARM200 and CERN COMPASS experiments which are looking for 
pentaquark and tetraquark states in pp collisions. Similar techniques can be used for searches 
at the RHIC collider. Charm exotics as well as hadrons with two or three quarks can be 
well produced at this energy. Hence, truly heavy ion collisions at RHIC opens the unique 
and tantalizing opportunity of producing and detecting multiply charm objects. They can 
possibly be detected by looking for a cascade of several weak decays of the charm quarks. 
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TABLES 



TABLE I. 


Hadron masses 


in the modified bag 


model (B 1 / 4 = 0.1684 GeV, a c 


= 0.4764, 


Zq = 1.585, m s 


= 0.342 GeV). 








TT\ £) T*t IP (3 

JJdl tide 


C[UclI ivo 


" 4 cxp L^evj 


''Hh L" fc!V J 


/\ 771 1 TO 

L\Ulj //icxp 




qqq 


n qiq 


D Q1Q 


fit 

lit 


p, u 


qq 


n 77fi 


D 77fi 


fit 
lit 


A 


qqq 


1 919 


1 919 


fit 
11 L 


7T, r/ 


qq 




n ^99 


u.u/o • 


o- 


coc 

odd 


X .u / o 


1 fi7^ 

1 .U ( o 


fit 

lit 


K 


sq 




U.OUO 


-t-lZr.^/O 


K* 


sq 


0.894 


0.920 


+2.9% 


A 


qqs 


1.116 


1.111 


-0.4% 


£ 


qqs 


1.191 


1.153 


-3.2% 


s* 


qqs 


1.384 


1.381 


-0.2% 


CI 


qss 


1.318 


1.302 


-1.3% 


"* 


qss 


1.533 


1.528 


-0.3% 




ss 


1.019 


1.060 


+4.0% 




TABLE II. 


Hadron masses with charm (m c = 1.788 GeV). 




JJdl tide 


nna v* Ire 


m exp [GeV] 


ma, [GeV] 


/\ 771 1 T71 

LAfU/ A/^cxp 


A 


/^i /^i r* 

qqc 


2.29 


2.29 


fit 

11 L 




cq 


1.87 


1.82 


-9 8% 




cq 


2.01 


2.01 


U. 1 /o 


7> s 


cs 


1.97 


1.98 


4-fl 7% 
TU. I /o 


D* 
u s 




2.11 


2.14 


-1-1 3% 


Vc 


CC 


2.98 


3.05 


4-9 4% 


T /?/> 
J ftp 


CC 


3.10 


3.15 


4-1 Q% 

t J- -y /o 


y, 


qqc 


2.45 


2.42 


-1 9% 




qsc 


2.47 


2.48 


+0.3% 


n c 


ssc 


2.70 


2.73 


+0.9% 


^cc 


qcc 


? 


3.66 






sec 


? 


3.82 






ccc 


? 


4.98 




s* 


qqc 


2.53(?) 


2.53 


-0.1% 


"* 
'— 'c 


qsc 


2.64 


2.67 


+1.0% 
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TABLE III. Hadron masses with beauty (m& = 5.30 GeV). 



particle 


quarks 


m exp [GeV] 


m th [GeV] 


Am/m cxp 


T(l-) 


bb 


9.46 


9.33 


-1.4% 


B 


bq 


5.28 


5.21 


-1.2% 


B* 


bq 


5.32 


5.31 


-0.2% 


B s 


bs 


5.37 


5.37 


fit 


b; 


bs 


5.42 


5.45 


+0.5% 




qqb 


5.64 


5.67 


+0.6% 


Vb 


bb 


? 


9.26 





TABLE IV. Charmed candidates and their main color magnetic and color electric interactions 
(q stands for u,d, and s quarks). 



Candidate 


V M /M 


Vk/Egc 


V E /E CC 


decays to 


V M /M 


vyE qc 


V E /E CC 


c5q [3] 


-16 


-16/3 





cqq + qqq 


-16 


-16/3 





cc4q [3] 


-16 


-16/3 


-8/3 


ccq + qqq 


-8 


-16/3 


-8/3 


cc4q [6] 


-28/3 


-40/3 


+4/3 


ccq + qqq 


-8 


-16/3 


-8/3 


3c3q [1] 


-8 





-8 


ccc + qqq 


-8 





-8 


3c3q [8] 


-14 


-12 


-2 


ccc + qqq 


-8 





-8 


3c6q [1] 


-24 





-8 


ccc + 2qqq 


-16 





-8 


3c6q [8] 


-46/3 


-12 


-2 


ccc + 2qqq 


-16 





-8 


4c2q [3] 


-8 


-16/3 


-8 


ccc + cqq 


-8 


-16/3 


-8 


4c2q [6] 


-4/3 


-40/3 


-4 


ccc + cqq 


-8 


-16/3 


-8 


4c5q [3] 


-16 


-16/3 


-8 


ccc+cqq+qqq 


-16 


-16/3 


-8 


4c5q [6] 


-12 


-40/3 


-4 


ccc+cqq+qqq 


-16 


-16/3 


-8 


5cq [3] 





-16/3 


-32/3 


ccc+ccq+qqq 





-16/3 


-32/3 


5c4q [3] 


-16 


-16/3 


-32/3 


ccc+ccq+qqq 


-8 


-16/3 


-32/3 


6c3q [1] 


-8 





-16 


2ccc + qqq 


-8 





-16 


6c6q [1] 


-24 





-16 


2ccc + 2qqq 


-16 





-16 


6c [11 








-16 


2ccc 








-16 
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TABLE V. Anticharmed candidates and their main color magnetic and color electric interac- 
tions (q stands for u,d, and s quarks). 



Candidate 


Vm/M 


vyE qc 


Ve/E cc 


decays to 


Vm/M 


VElE qc 


V E /E C 


c4q [3] 


-16 


-16/3 





D + qqq 


-8 


-16/3 





ccqq [3] 


-8 


-16/3 


-8/3 


2D 





-32/3 





ccqq [6] 


-4/3 


-40/3 


+4/3 


2D 





-32/3 





cc5q [3] 


-16 


-16/3 


-8/3 


2D + qqq 


-8 


-32/3 





cc5q [6] 


-12 


-40/3 


+4/3 


2D + qqq 


-8 


-32/3 





ccc3q [1] 


-8 





-8 


ccc + qqq 


-8 





-8 


ccc3q [8] 


-14 


-12 


-2 


ccc + qqq 


-8 





-8 


ccc6q [1] 


-24 





-8 


ccc + 2qqq 


-16 





-8 


ccc6q [8] 


-46/3 


-12 


-2 


ccc + 2qqq 


-16 





-8 


4cq [3] 





-16/3 


-8 


D + ccc 





-16/3 


-8 


4c4q [3] 


-16 


-16/3 


-8 


D+ccc+qqq 


-8 


-16/3 


-8 


4c4q [6] 


-28/3 


-40/3 


-4 


D+ccc+qqq 


-8 


-16/3 


-8 


5cqq [3] 


-8 


-16/3 


-32/3 


2D + ccc 





-32/3 


-8 


5c5q [3] 


-16 


-16/3 


-32/3 


2D+ccc+qqq 


-8 


-32/3 


-8 


6c3q [1] 


-8 





-16 


2ccc + qqq 


-8 





-16 


6c6q [1] 


-24 





-16 


2ccc + 2qqq 


-16 





-16 
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TABLE VI. Estimate of the binding energy of charmed candidates assuming M = 18 MeV, 
E qc = 50 MeV, E cc = 70 MeV (q stands for u,d, and s quarks). 



Candidate 


AVm/M 




AV E /E CC 


Energy gain [MeV] 


cc4q [3] 


-8 








144 


cc4q [6] 


-4/3 


-8 


+4 


144 


3c3q [8] 


-6 


-12 


+6 


288 


3c6q [1] 


-8 








144 


3c6q [8] 


+2/3 


-12 


+6 


168 


4c2q [6] 


+20/3 


-8 


+4 





4c5q [6] 


+4 


-8 


+4 


48 


5c4q [3] 


-8 








144 


c4q [3] 


-8 








144 


ccqq [3] 


-8 


+16/3 


-8/3 


64 


ccqq [6] 


-4/3 


-8/3 


+4/3 


117 


cc5q [3] 


-8 


+16/3 


-8/3 


64 


cc5q [6] 


-4 


-8/3 


+4/3 


165 


ccc3q [8] 


-6 


-12 


+6 


288 


ccc6q [1] 


-8 








144 


ccc6q [8] 


+2/3 


-12 


+6 


168 


4c4q [3] 


-8 








144 


4c4q [6] 


-4/3 


-8 


+4 


144 


5cqq [3] 


-8 


+16/3 


-8/3 


64 


5c5q [3] 


-8 


+16/3 


-8/3 


64 


6c6q [1] 


-8 








144 
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TABLE VII. Charmed candidates and their masses in the modified bag model (q denotes the 
up or down quark, a star indicates that the threshold mass is known). 



candidate 


charge 


mass [GeV] 


estimated thresl 


csudud 


+1 


3.58 


3.40* 


cssudq 


0,+l 


3.66 


3.58* 


ccsudq 


+2,+l 


4.87 


4.76 


ccssud 


+1 


5.04 


4.94 


cccsud 


+2 


6.29 


6.10 


cccssudud 


+2 


7.35 


7.21 


4cud 


+3 


7.53 


7.27 


4csq 


+3,+2 


7.72 


7.45 


4csudud 


+3 


8.70 


8.38 


4cssudq 


+3,+2 


8.84 


8.56 


5cq 


+4,+3 


9.00 


8.64 


5cs 


+3 


9.14 


8.80 


5csudq 


+4,+3 


10.0 


9.76 


5cssud 


+3 


10.2 


9.92 


6c 


+4 


10.3 


9.96 


6csud 


+4 


11.5 


11.1 


6cssudud 


+4 


12.6 


12.2 



[GeV] 



TABLE VIII. Anticharmed candidates and their masses in the modified bag model (q denotes 
the up or down quark). 



candidate 


charge 


mass [GeV] 


threshold 


csudq 


o,-i 


3.13 


2.91 


cssud 


-1 


3.31 


3.09 


ccud 


-1 


4.04 


3.73 


ccsq 


-1-2 


4.23 


3.84 


ccsudud 


-1 


5.24 


4.78 


ccssudq 


-1,-2 


5.38 


4.88 


4cq 


-2,-3 


7.25 


6.85 


4cs 


-3 


7.39 


6.95 


4csudq 


-2,-3 


8.33 


7.57 


4cssud 


-3 


8.51 


7.67 


5cud 


-3 


9.27 


8.71 


5csq 


-3,-4 


9.49 


8.82 


5csudud 


-3 


10.4 


9.44 


5cssudq 


-3,-4 


10.6 


9.54 



[GeV] 
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FIGURES 
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fc 

FIG. 1. Binding energy of charmed matter for various strong coupling constants. For compar- 
ison the case for strange matter indicated by m s is also shown. 
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FIG. 2. Allowed region for the formation of charmed matter as a function of strange and 
charmed quark fraction, f s and f c respectively. The different shaded regions correspond to different 
values of the strong coupling constant a c . The region labeled 'fit' corresponds to our bag fit of 
a c = 0.476 and B 1 / 4 = 168 MeV. 
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FIG. 3. Production probabilities for multiply charmed quark exotics at RHIC (s 1 / 2 = 200 GeV) 
as a function of the coalescence parameter po with shadowing and a quench. 
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FIG. 4. Production probabilities for multiply charmed quark exotics at RHIC as a function of 
the coalescence parameter po for the unquenched case with no shadowing. 
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FIG. 5. Probability density of a quark of flavor q(u, d), s, c, b for a bag of radius R = 1 fm as a 
function of the scaled radius. 
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FIG. 6. Electric interaction Ey between two quarks with flavor i, j = q(u, d), s,c as a function 
of the bag radius R. 
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FIG. 7. Magnetic interaction My between two quarks with flavor i,j = q(u, d), s,c as a function 
of the bag radius R. 
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